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Abstract 

Within the i^MSSM, a Minimal Supersymmetric neutrino See-saw Model, Lepton 
Flavour Violating Higgs couplings are strongly enhanced at large tan /3 ( > 30) , 
which can lead to BK{H^/A^ Tfi) ~ 0(10-'^), for Mj^o/^o > 160 GeV. Enhance- 
ments on the production of Higgs bosons, through the gluon fusion mechanism, 
gg H^/A^, and the associated production channel gg,qq — > bbH^/A^, whose 
rates grow with tan /3, as well as the mass degeneracy that occurs between the 
and A^ states in this regime, also contribute to further the possibilities to detect a 
heavy Higgs signal into Tfi pairs. We show that the separation of Tfi Higgs events 
from the background at the upcoming CERN Large Hadron Collider could be done 
for Higgs masses up to about 600 GeV for 300 fb~^ of luminosity, for large tan/? 
values. However, even with as little as 10 fb~^ one can probe /A^ masses up to 
400 GeV or so, if tan/3 = 60. Altogether, these processes then provide a new Higgs 
discovery mode as well as an independent test of flavour physics. 



1 Introduction 



It is possible that some of the most exciting times in particle physics will come soon, as 
the upcoming Large Hadron Collider (LHC) will enable us to probe the mechanism of 
Electro- Weak Symmetry Breaking (EWSB). The standard picture contains a single Higgs 
state that couples to the other fundamental particles with an intensity proportional to 
their masses [1]. In the case of fermions, the Standard Model (SM) Higgs couplings are 
diagonal in flavour space, due to the fact that the Higgs couplings and the fermion mass 
matrices are both diagonalised by the same bi-unitary rotations. However, this picture 
ceases to remain valid in many extensions of the SM. For instance, in the general Two 
Higgs Doublet Model of Type III (THDM-III), where both Higgs doublets couple to 
both types of up- and down-type fermions, there appear non-diagonal Higgs couplings, 
which lead to interesting Lepton Flavour Violating/Flavour Changing Neutral Current 
(LFV/FCNC) Higgs phenomenology [2]. In turn, even though the Minimal Supersym- 
metric Standard Model (MSSM) is a Type II Two-Higgs Doublet Model (THDM-II) 
at tree level - with additional mass and coupling relations enforced by Supersymmetry 
(SUSY) - this structure is not protected by any symmetry, so that loop effects can effec- 
tively render it a THDM-III. In addition, the detection of neutrino oscillations [3, 4] seem 
to suggest that there is a large mixing between the second and third families in the lepton 
sector, which could also appear in new scenarios that are contained in some extensions of 
the SM [5]. In particular, within SUSY models, the pattern of LFV effects at the Planck 
or Grand Unification Theory (GUT) scales could be reflected in the structure of the soft 
SUSY-breaking terms, i.e., in the slepton mass matrices, which in turn can communicate 
these to the Higgs sector through radiative effects [6] . 

Detectable effects of LFV Higgs couplings could show up in the decay r 3/^, which 
is a particularly sensitive probe at large tan/? [7] (the ratio of the two Higgs vacuum 
expectation values in the MSSM), with a Branching Ratio (BR) scaling as tan^/?, a 
phenomenon which may render this mode detectable at the LHC. The relevance of the 
LFV Higgs decay cj) — > t/j,^ for Higgs phenomenology at the LHC was discussed in Ref. 
[8]^, within the context of several extensions of the SM. In particular, it was shown that 
a large BR(0 r^) (of order 0.001-0.01) could easily be achieved in the THDM-III. 
Moreover, it was shown there that large LFV Higgs couplings were not in conflict with 
any low energy constraints, such as LFV decays of r's. Calculations of the SM Higgs 
BR(0 — > T//) showed it to be be very suppressed (< 10"-*^^) whilst in the (constrained) 
MSSM the corresponding rates could be enhanced for some of the Higgs states. 

In this paper we are interested in discussing further aspects of LFV phenomenology 
entering the Higgs sector of the MSSM, by investigating the possibility of a new detection 
mode for heavy H^/A^ Higgs bosons {Mho/a" ^ 2M^y± GeV) at the LHC. In particular, 
we shall demonstrate that the decays of the heavy neutral Higgs bosons of the MSSM to 
T/i pairs are sizable and represent a very sensitive probe of LFV physics. We calculate 

^Hereafter, unless otherwise specified, the label ^ will refer to a generic Higgs state. 
^However, a calculation of the actual LFV Higgs decay rates was presented first in [9] . 
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the rates for the BR(/J°//1° Tfi), involving the heavy neutral MSSM CP-even Higgs 
state and the CP-odd one, and find that they can be as large as 10""^, for high values 
of tan/3, while those for the BR(/i° — > r//), involving the light neutral MSSM CP-even 
Higgs state, remain rather small in comparison. Furthermore, it should be recalled that 
MSSM Higgs production at hadron colliders is enhanced for large tan P, both via gluon 
fusion and in association with bb pairs. Moreover, in this tan f3 regime, there appears a 
degeneracy for the masses of the and states, which essentially doubles the event 
rate of the overall Higgs signal. All such effects enable one then to reach detectable 
levels at the LHC for H^/A^ — > t/j, signals. Finally, by extending previous studies on 
the signal-to-background separation at the LHC, we show that detection of H^/A^ LFV 
Higgs decays into tau-muon pairs could be achieved for Higgs masses as high as 600 GeV. 
In short, these LFV decay modes can provide a new detection channel for the heavy Higgs 
bosons of the MSSM, which would in turn give not only important evidence for SUSY 
but also, along with the modes S° — > /i/i, r — > S/i, t ^ /I'j and /i — > 67, probe the form 
of the neutrino Yukawa mass matrix. 

The plan of the paper is as follows. The underlying aspects of the model and the 
calculation are laid out in Sect. 2. The LFV Higgs signals are characterised in Sect. 3. 
The numerical calculation of the cross sections and decays for the relevant Higgs modes is 
pursued in Sect. 4, including the determination of signal-to-background event rates and 
the proof of dctcctabihty of LFV Higgs signals at the LHC. Summary and conclusions 
are found in Sect. 5. 

2 Slepton Mixing and LFV in the Higgs Sector 

One of the most attractive explanations for the observed neutrino masses [3, 4] is the 
"see-saw" mechanism [10], which includes Dirac masses {uid) as well as Majorana masses 
(Mr). Atmospheric neutrino data favours a Ur mass of about 0.04 eV [11]. Thus, for 
Dirac neutrino masses of the order of the corresponding up-quark masses, i.e. ~ 
100 — 200 GcV, as predicted in a GUT such as SO(IO), one finds that the right-handed 
Majorana mass, Mji, needs to be of order 10^^ GeV. Majorana neutrino masses imply 
LFV within the Minimal Supersymmetric see-saw Model, which is defined as the Minimal 
Supersymmetric Standard Model augmented by three heavy right-handed neutrinos, i/r^, 
LFV interactions can be communicated directly from z//j's to the sleptons and from these 
to the charged leptons and Higgs bosons. The initial communication takes place through 
renormalisation group fiow of the slepton mass matrices at energies between Mpianck and 
Mr. The presence of ur states at scales above Mr leaves an imprint on the mass matrices 
of the sleptons, which propagates down to the Electro- Weak (EW) scale. This effect has 
been used to predict large BRs for r — > /i7 and — ^ 67 within the MSSM [12, 13, 14]. 
To derive the effective Lagrangian for the LFV lepton-Higgs interactions, we begin 

^Henceforth, we will use the notation z/MSSM to indicate such an extension of the MSSM. 
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with the Yukawa Lagrangian: 



-C = iRYiLLHd + VrY^LlH^ + \vIMr vr, (1) 

where Ir, Ll and vr represent the right-handed charged leptons, left-handed lepton 
doublets and right-handed neutrinos, respectively, while denote the Higgs doublets 

of the MSSM. YJ, and Mr are 3x3 matrices in flavour space. We shall work in a basis 
in which both Yi and Mr have been diagonaliscd, but where Yi, remains an arbitrary 
complex matrix. Lepton number is violated in this Lagrangian due to the presence of 
the vr Majorana mass term. 

Furthermore, the 6x6 slepton mass matrix is written in terms of 3 x 3 blocks, as 
follows: 



Ml 



ML ML 



(2) 



where 



ML = M| + Mf + M|o cos 2/5 (Tg^ -Q^sin^^vi/) (3) 
MIr = Ml + Mf + Mlo cos 2PQ I sin^Ow (4) 
ML = Aiv cos /3/ V2- Ml lutein 13 (5) 

Here v = 246 GeV, 9w is the weak/Weinberg angle while and Ai denote the higgsino 
mass parameter and trilinear slepton couplings, respectively. Myy±^zo are the masses of 
the W^, Z° gauge bosons and Mi the lepton mass matrix. We will work in a basis where 
Ml is diagonal. 

When the SUSY-breaking slepton mass matrix {M^ij evolves from the scale M at 
which flavour-blind SUSY-breaking is communicated to the visible sector, down to the 

slepton mass scale (assuming M > Mr), one obtains a flavour-mixing piece that 
corrects the slepton soft mass terms, i.e., M| — > M| -|- AM|, with the latter given by: 

where mo is a common scalar mass evaluated at the scale Q — M and i ^ j- If one 
assumes that the A-terms are proportional to Yukawa matrices, then: 

(AM|),,=.£(yJn),,. (7) 

where 

(8) 

and a is 0{1). In the simplest SUSY-breaking scenarios, gravity plays the role of messen- 
ger and M — Mpianck- Global fits to neutrino data favour large mixing between i/^ and 
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Vt- and also between u^. and [11]. Thus, we shall consider here a form for niy with 0{V) 
entries in the 23 and 32 elements [15]. If we further assume that {Mji)ij is an identity 
matrix, then {YjYi,)23 will also be of 0(1). 

This source of LFV interactions can be transmitted to the Higgs sector as well, because 
radiative effects could induce flavour mixing in the Higgs couplings. These corrections 
allow the neutral Higgs bosons to mediate FCNCs, in particular /i/i [6] can reach 

BRs at large tan/3 that can be probed by Run II of the Tevatron [16]. For the leptonic 
sector, the Feynman graphs that induce such corrections involve loops of sleptons and 
charginos/neutralinos, which are transmitted further to induce LFV Higgs couplings. 
To derive the SUSY-induced THDM-III, one can write an effective Lagrangian for the 
couplings of the charged leptons to the neutral Higgs fields, namely: 

-C = IrYiElH^, + IrYi (eil + e^Y^) ElH^^* + h.c. (9) 



where ei,2 include the slepton radiative effects. LFV couplings results from our inability 
to simultaneously diagonalise the term Yi and the non-holomorphic loop corrections, 
e2YiYjY^. The contributions from higgsinos and gauginos, which are approximated as 
mass eigenstates, can be written as follows, 

~ ^c/^M/i (/^', ^L' ^L' ^i) ' (10) 

in which £a = jl, e, i>i, and — f, Or- Mj = Mi 2 are the U(l) and SU(2) gaugino masses, 
while ^ is defined in eq. (8). For our purposes, the function /j can be evaluated in the 
limit a = b = c = d, for which fi{a, a,a,a) = 1/ (6a^). 

Since the charged lepton masses cannot be diagonalised in the same basis as the 
Higgs couplings, this will allow neutral Higgs bosons to mediate LFV processes, with 
rates proportional to e|. The term proportional to ei will generate a mass shift for the 
charged leptons that will appear as a second-order effect [17]. 

Then, the LFV interactions relevant for the Higgs sector phenomenology can be writ- 
ten in terms of Higgs mass eigenstates as follows (0^ = h^, H^, A^): 

where 



A# ^ ... (12) 



-e2tan/3p^ (yJY,), 
sin f3r]^ 

and {Tih,riH,'qA) = (sin a, — cos a, sin {ph, Ph, Pa) = (cos(/9— a), — sin(/5— a), — i), with 
a denoting the Higgs mixing angle. Constraints on the LFV (r^j /i^) -Higgs interaction 
can be obtained from the LFV r-dccays (e.g., r 3fi), which can be generated via 
exchange of h^, and . For instance, for the case in which p = Mi = M2 = = rrij), 
Mr = 10^^ GeV and {YJY^)s2 = 1, Ref. [7] finds that £2 ~ 4 x 10-^ which is stable with 
respect to changes in the SUSY spectrum. Then BR(t 3p) ~ (1 x 10"'^) x (tan P/60f x 
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Figure 1: The variation of the LFV BR of heavy Higgs bosons function of 

Mj^o{Mho = M^o) for for tan/3 = 30,40 (with set A) and 50,60 (with set B). 

(100 GeV/M^o)^, which puts the r — > mode into a regime that is experimentally 
accessible at current B-factories. At the LHC and SuperKEKB, hmits in the region of 
10~^ should be achievable [18], allowing an even deeper probe into the model parameter 
space. On the other hand, Tevatron has already constrained the large tan/5 domain, 
from the search for the decay Bg fifi [19], that seems to exclude the value tan/3 = 60, 
which is preferred by the requirement of Yukawa unification. However, it is possible to 
evade such constraints, for instance in SUSY breaking scenarios where slepton and squark 
masses do not have any strong correlation. In any case, given that such limits depend 
on multiple MSSM parameters, which makes it difficult to draw a general conclusion, it 
is certainly preferable to test the resulting LFV Higgs couplings directly at the LHC, as 
it is discussed in the next section. 



LFV Higgs decays have been evaluated within the general MSSM with a particular ansatz 
for the trilinear A terms in [5], where was found that a BR(^° — > r/x) ~ 10~^ — 10"''' 
could well be achieved. Subsequently, Refs. [20, 21] and [22] presented a more detailed 
calculation of the BK{h^ rfi), within the MSSM, which was essentially in agreement 
with the previous result. In fact, Ref. [22] also reported a complete one-loop calculation 
of the LFV Higgs decay within the SM extension with massive neutrinos, using a realistic 
pattern of neutrino masses and mixings, which resulted in very suppressed LFV Higgs 
decays (with BR of order 10~^° or less). Afterwards, [23] presented a detailed study of 
the prospects to detect LFV Higgs decays at the Tevatron and LHC, concluding that it 



3 The LFV Higgs Decays H^A 
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is certainly possible to detect such decays in the THDM-III within the Higgs mass range 
114-160 GeV, approximately. Later on, Ref. [24] presented a more realistic study of the 
signal and backgrounds, essentially reaching the same conclusions. Subsequently, it was 
also studied in detail the mass-matrix ansatz used in the THDM-III, by Ref. [25] , while 
the evaluation of the corresponding LFV Higgs decays was presented in [26]. Mixing of 
the SM fermions with other exotic fermions was also shown to be a possible source of 
LFV in the Higgs sector [27], resulting in BRs of the order 0.01-0.001 again, which could 
clearly be detectable too. Bounds on LFV Higgs decays at the Tevatron were reported 
by the CDF collaboration in [28]. 

In this paper we shall concentrate on the LFV heavy Higgs decays to r^u, which has a 
very small BR within the context of the SM with light neutrinos, so that this channel is 
potentially an excellent window for probing new physics. Although the previously men- 
tioned works have studied LFV Higgs decays within the i/MSSM, the specific evaluation 
of the BRs of the LFV decays for /A^ in the heavy mass range has not been studied. 
Hence a discussion of the corresponding detectability at the LHC has not been presented 
so far either. To remedy this is our aim in the present paper. 

In order to derive the formulae for the LFV Higgs decay widths, we notice that the 
quantity inside the square brackets in eq. (11) corresponds to the Higgs-lepton coupling 
0feTT, which will be denoted by gtt>TT- Thus, we can write the LFV Higgs coupling 0fcT// 
as 

g4>Tti = g<t>TT^tj- (13) 

The decay width for the generic process </)fc (in which wc add both final states 

r^/i^ and t~ij,~^) can then be written in terms of the Higgs decay width r(0 —>■ tt), as 
follows: 

r(0^r/x) = 2\\tfT{^^TT), (14) 

so that the LFV Higgs BR can in turn be approximated by BR(0 — > r/x) = 2|A^^pBR(0 — > 
tt). 

We are interested in studying the large tan/3 domain (i.e., (3 7r/2), where the LFV 
Higgs couplings are enhanced. It is also simpler to work in such so-called decoupling 
regime of the MSSM Higgs sector, which in fact is quite general since it is reached even for 
moderate values of M^o (— 200 GeV). In this case we have that A^^ — >• 0, €2 tan /3, €2 tan (3 
for (j)k — hP,H^,A^, respectively. Therefore, the LFV decays of the hght Higgs boson 
{hP) are suppressed for most regions of parameter space. Conversely, for the above 
mentioned choices of SUSY parameters yielding 62 = 4 x 10""^ (which we call set A) 
and with MuQ j^o 160 GeV, one obtains BR(i/°/A° tt) ~ 0.12, which in turn gives 
BR(if'^/y4'^ — > r/i.) ~ 2.9 x 10~^ for tan j3 = 30. For another set of parameters with a large 
fjL hmit, i.e. /j, » Mi^2 (which we call set B), one gets 62 ~ 8 x 10~^, which will produce 
a larger BR for H'^ /A^ — > r/x. More in general, we calculate the LFV Higgs decay rates 
in the channels H^/A^ ^ r^u - by appropriately modifying the HDECAY program [29] and 
using the formula in eq. (14) - as a function of the Higgs masses and corresponding LFV 
Higgs couplings. 
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Figure 2: The dependence of number of events of the LFV Higgs signals at the LHC on the 
degenerate Higgs masses M^jo and M^o, including the sum of production cross sections cr{gg — > 
H^) + cr{gg — > A^) times the corresponding decay rates BR(i7'^ — > r/i) and BR(A'^ — > r/i), for 
two representative values of tan/3 = 30 (+set A) and 60 (+set B). We are assuming a detection 
efficiency of 3% and two values of integrated luminosities 30 fb~^ and 100 fb~^. 



It is appropriate to mention at this point that reference [24] did include a discussion 
of the LFV H^/A^ decays. However, the authors concentrated on the mass range below 
2M\Y± ^ 160 GeV. They hint in fact that above this mass range the modes H^/A^ t/j, 
will be suppressed because the channels W~^W~ and Z^Z^ would be open and dominate 
the total decay width. However, this is not true. The reason is twofold. Firstly, the 
A^ - being a CP-odd state - does not couple to vector boson pairs. Secondly, although 
the state does couple to W~^W~ and Z^Z^ pairs, in the heavy Higgs mass limit such 
coupling is considerably suppressed. Our calculation takes correct care of these aspects. 



4 Signal-to-Background Analysis 

Once one folds in the values of LFV Higgs BRs with the main heavy neutral Higgs 
production modes at the LHC, it is clear that LFV Higgs decays into r/i pairs may be 
detectable at the CERN hadron collider. As a benchmark, according to previous studies 
[24], with SM-like cross sections and < 160 GeV, one could detect at the LHC the 
aforementioned LFV Higgs decays with a BR of order 8 x 10~^. We intend to push forward 
the region of detectability into higher mass values within the MSSM, by exploiting the 
aforementioned tan (3 enhancement and the fact that our LFV Higgs BRs become almost 
constant for heavier Higgs masses. Besides, for heavier Higgs masses one should expect 
a much larger background reduction, as compared to the lower Higgs mass case, owning 
to the much harder energy spectra for the emerging r- and /i-leptons. 

As already mentioned, high values of tan /3 are also associated with a large 6-quark 
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Figure 3: The dependence of number of events of the LFV Higgs signals at the LHC on the 
degenerate Higgs masses Mfjo and M^o, including the sum of production cross sections a{pp — > 
bbH^) + (7{pp — > bbA^) times the corresponding decay rates BK{H^ Tfi) and BR(A'' — > Tfi). 
Other parameters are the same as in Figure 2. 



Yukawa coupling, which in turn can produce an enhancement of the Higgs production 
cross sections at hadron colliders via both gluon fusion and associated production with b- 
quark pairs: i.e., gg —>■ H^/A^ (via triangle loops at lowest order) and gg, qq — > bbH^/A^ 
(at tree-level), respectively. (We calculate these two production processes here by using 
the HIGLU and HQQ programs in default configurations [30].) Figure 1 shows the relevant 
LFV Higgs BRs as a function of M^o . The LFV BR is basically the same for both if ° and 

and its variation with M^o is very mild. In Figures 2 and 3 we present the expected 
LFV Higgs event rates as a function of degenerate Higgs mass M^o for two values of 
tan P = 30 and 60 and two values of LHC luminosities, L = 30 fb~ and 100 fb~ , 
respectively. It is clear from these two plots that LFV Higgs rates con be substantial 
even at large Higgs masses. 

In Ref. [23] it was proposed a series of cuts to reconstruct the hadronic and electronic 
r decays from (p —>■ Tfi and separate the signal from the background, which is dominated 
by r-pair production via Drell-Yan modes (i.e., qq 7*,^*^*^ — r+r~) and qq,gg 
W~^W~ T^UrT~Vr. In fact, it should be recalled that the decay product distributions 
of r-leptons generated in the decay of Higgs bosons are notably different from those 
emerging in gauge boson decays, because of the different spin of the primary objects. A 
more realistic search strategy for the LHC based on the cuts of Ref. [23] was presented 
in Ref. [24], where one can find detection efficiencies in the Higgs mass range 120-160 
GeV. The typical figure goes from about 2% for M/^o = 120 GeV up to about 3% for 
Mffi = 160 GeV, where it starts stabilising. Although one expects that this detection 
efficiency will increase for heavier Higgs masses, in order to use a conservative estimate, 
we shall use the 3% figure throughout in our estimates in the remainder of the paper. 
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Figure 4: The 95% CL exclusion and 5o" discovery reaches for the LFV Higgs signals at the 
LHC as a function of the degenerate Higgs masses Mj^o and M^o and the integrated luminosity, 
including the sum of production cross sections cr{gg H^) and (T{gg — > A9) plus cr{pp — > bhH^) 
and a{pp hhA^) times the corresponding decay rates BR(//'^ — > t^) and BR(A'^ — > Tfi), for 
tan/3 = 30 (+set A) (dashed lines) and 40 (+set A) (solid lines). 
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Figure 5: Same as Figure 4, but for tan/3 = 50 (+set B) (dotted-dashed hues) and 60 (+set B) 
(solid lines). 
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Adopting the background rates estimated in Ref. [24] for Higgs masses up to 200 
GeV and trivially extrapolating them to heavier Higgses, we are then in a position to 
compare the yield of our signals (see Figures 2-3) with that of the total background, 
thereby estimating both a 95% Confidence Level (CL) exclusion hmit and a 5a discovery 
reach. The scope of the LHC in both respects, as a function of the Higgs masses and 
machine luminosity, is then well described by Figures 4 and 5. To display our results, we 
have choosen the combinations tan /? = 30 and 40 with SUSY parameters of set A, and 
tan 13 — 50 and 60 with SUSY parameters of set B. (Other combinations should lay within 
these results.) Clearly, with an integrated luminosity of 300 fb~^, for tan/3 = 30, one 
can detect a signal for Higgs masses up to about 260 GeV, which is already significantly 
above the 2M^± mark of Ref. [24], while, for tan (3 = 60, Higgs masses up to even 600 
GeV can be probed. However, since for tan (3 = 60 one is dangerously close to the bounds 
from B ^ nn, the reader may well refer instead to those for tan P — AO (with SUSY 
parameters of set A) and tan/3 = 50 (with SUSY parameters of set B). In these cases 
too we find that it will be possible to extract a LFV Higgs signal at the LHC for heavy 
Higgs masses, up to about 300(500) GeV for tan/3 = 40(50) with 300 fb~^ of luminosity. 

5 Summary and Conclusions 

We have demonstrated that LFV effects in the slepton sector of the i/MSSM can generate 
LFV couplings between leptons and neutral Higgs bosons leading to large BRs for LFV 
Higgs decays into lepton pairs. In particular, we have calculated the BRs of the processes 
H^/A^ TyU and found that they can be as large as 3 x 10"*^, while the BR(/i° Tfi) 
is only about ~ 10^*^. Furthermore, these rates occur for large values of Mf^o/j^o and 
tan /3, a configuration also responsible for a strong degeneracy between the masses (and 
couphngs) of the and ^4° states, producing an overall Higgs event rate which is double 
the one of cither Higgs state alone. These LFV Higgs modes can be extracted at the 
LHC for Higgs masses slightly beyond 600 GeV, provided tan/3 = 60. For smaller values 
of this parameter though the LHC scope greatly diminishes, reducing to just above 260 
GeV in Higgs mass for tan/3 = 30. These values can only be reached at 300 fb~^ of 
luminosity. However, even with a modest 10 fb~^, one could probe Higgs masses up to 
400(415) GeV, provided tan/3 = 50(60). Besides, in view of the assumptions made on 
detection efficiencies for the signal, we believe these conclusions to be rather conservative. 
Altogether, these novel channels complement the modes fifi, r 3/i, r ^ fi'-f and 

— > 67 in order to provide evidence for SUSY and key insights into the form of the 
neutrino Yukawa mass matrix. More detailed experimental simulations are now awaited. 
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